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Abstract 

We outline the basic ideas of the topologigical mechanisms of su- 
perconductivity. A gauged model of correlated electronic system where 
a topological fluid is formed as a result of a strong interaction is dis- 
cussed. 



1 



1 Citerion for Superconductivity 



1.1 London Equations as Hydrodynamics of an Ideal 
Liquid 

Our understanding of superconductivity starts from the London equations. 
They describe the hnear response of superconductive matter to an external 
electromagnetic field: 

Wxj = -^B^^\ (1) 

Vp + vo-'dtj = KE''''' . (2) 

The first equation implies the Meissner effect (A is the London penetration 
depth). The second means that the matter is compressible Q (k is a com- 
pressibility) and implies an ideal conductivity {(j{uj) ~ iuj~^). Of course, due 
to gauge invariance (or the Maxwell condition V xE = dtB), the penetration 
depth and compressibility are related by k = (Afo)~^. 

These two equations may be written as one. In the Lorentz gauge VA + 
v^'^dtAo = 0, for example, one has 

J = -A-'l'="* . (3) 

These equations describe the hydrodynamics of the superconductor at 
low temperatures. 

If we replace E*^^^ in (|^) by the gradient of pressure Vp, and i?^^* in (|^) 
by a torque, we obtain the linearized Euler equation for an ideal compressible 
liquid. Therefore the Meissner effect is a direct consequence of the fact that 
superconductive matter is a liquid. 

Let us elaborate this simple but important observation a bit further. The 
London equations can also be written in terms of current-current correlators 
at small u and k. Varying these equations over A^^* and taking into account 
the continuity equation dtp + Vj = 0, we find 

{p{k, uj)p{-k, u)) = -K— 2— (4) 

— v^k-^ 

^ As a result of Coulomb interactions the liquid is in fact incompressible. Below we 
neglect the Coulomb interaction between charge carriers and consider the electromagnetic 
field as external. 
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and 

{j±{k,uj)j^{k,uj)) = Vjl = , (5) 

where j± is the transverse component of the current. Now it is obvious that 
our matter is a hquid. The second equation (which equivalently reflects the 
Meissner effect) imphes that there is no gapless transverse mode (transverse 
sound), i.e. that there is no shear modulus - this is the definition of a liquid. 
A transverse sound is a feature of a solid for which, instead of (|), one has 

1 UJ^ 

{j±{k,uj)j^{k,u)) = -—-^—^ . (6) 

The density- density correlator (^ shows a gapless mode (longitudinal sound) 
which means that the liquid is ideal and compressible, i.e. may flow without 
dissipation. 

All these equations can be linked to a Hamiltonian: 

H = g[f + vlip - p)'] + A^^'j + Arip - P) (7) 

where p is an average density. This is the Hamiltonian of linear hydrody- 
namics. 



1.2 Landau Criterion 

To be a superfluid a system must have no dissipation.lt seems difficult to 
combine a soft mode in density modulation and the absence of dissipation. 
A Fermi liquid, for example, shows Landau damping, 

(i,ft.),,(*,.)> = -^^jJ— ^, (8) 

due to gapless particle-hole excitations {x and 7 are constants). To have no 
dissipation a liquid must have no gapless density modes rather than longitu- 
dinal sound. In particular the single particle spectrum of the homogeneous 
liquid must have a gap. In known examples it is sufficient that all other 
excitations, which change a local density also have a gap. This is the subject 
of the Landau criterion of superfluidity. In its grotesque form it states 
that 
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the spectrum must contain longitudinal sound and the single par- 
ticle spectrum must have a gap. 

It is believed that this condition on the spectrum of quantum hquid is 
sufficient to derive the Meissner effect. 

All of this is true for a homogeneous liquid. The true check whether a 
liquid is superconductive may be made only if some non-zero concentration 
of impurities does not result in resistivity. A further belief is that weak 
impurities do not lead to resistivity whenever the Meissner effect holds in 
the absence of impurities. 

The Landau criterion seems to be sufficient in spatial dimensions greater 
than one. In one dimension it fails for a very simple reason, namely, in one 
dimension we cannot distinguish between a liquid and a solid since there are 
no properties such as transversal modes and shear modulus. There is no 
Meissner effect either. As a result a single impurity pins down the flow in 
the same way as a single impurity pins down a solid. 

In the early days of superconductivity, Frohlich p| noticed that in a one- 
dimensional metal an incommensurate charge density wave (CDW) will slide 
through the lattice unattenuated. Since it carries an electric charge and 
since a gap has developed in the electronic spectrum, Frohlich concluded 
that the ground state of his system is superconductive. Although, a sliding 
charge density wave indeed contributes to conductivity, the Frohlich super- 
conductivity in one dimension was considered nothing more than a theoretical 
curiosity, because of a variety of pinning mechanisms 0. 

The failure of the charge density wave mechanism in one dimension does 
not devaluate Frohlich's ideas, which may be valid in higher dimensions [|]. 
The distinctive property of Frohlich superconductivity is a topological order - 
a topological character of the ground state. We therefore call it the topological 
mechanism of superconductivity. One of the best known examples of this 
phenomenon in two dimensions is the superconductivity of particles with 
fractional statistics or anyon superconductivity (see e.g.[0). 

There are two reasons why the topological mechanism of superconduc- 
tivity is of interest. First of all it is a fundamentally new mechanism of 
superconductivity and superfluidity. Secondly it has been found in models of 
strongly correlated electronic systems, namely in the doped Mott insulator 
(some people consider it to be a mechanism of high temperature supercon- 
ductivity) . 
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Below we discuss some models of topological superconductivity and dis- 
cuss why they are relevant to the physics of correlated electronic systems. 



1.3 Topological order 

Jumping ahead, let us discuss some features of the ground state with a topo- 
logical order. Let the wave function of the system be ip{{xi}) where {xi} are 
coordinates of particles. Its defines a fiber bundle. We shall refer the Chern 
number of this fiber bundle ClS db topological order. 

Consider for example a two dimensional system. A topological order 
implies that the ground state wave function has zeros with a non-vanishing 
signature. In this case the Chern's number is 



where the integral goes over the space. This equation can be also applied 
for 3-dimensional systems. In this case the integral goes over the boundary 
of the system. In one dimension, let us define the Chern number as follows. 
Consider a closed one- dimensional space as a boundary of a domain of the 
complex plane and the wave function as the value of some analytical function 
in the domain. Then (j^) gives the Chern number. 

Perhaps the simplest example were a topological order occurs is particles 
in a magnetic field. Then the Chern number is merely Hall conductance. 
However a topological order may occur in correlated electronic systems even 
without magnetic field, but as a result of a strong interaction. Thse sytems 
are not exotic. They were known from the early days of quantum physics in 
one dimension. They have also recently been discovered in quantum mag- 
netism 10, H, ^, [ll| and intensively discussed with respect to fractional 
statistics in two dimensions. A topological order also occurs in some models 
of electronic systems with strong interaction in two and three dimensions. 
Some people even think that a topological order is a universal feature of 
correlated quantum systems. 

An illustrative example is given by the simplest solvable model of particles 
in one dimension with a l/r^ interaction: 






(10) 
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The exact wave function of this model is 



i!{xi....XN) = Y[{Xi - Xj) 



A 



(11) 



This wave function tells us that each collision contributes vr to the phase shift. 
The Chern class, therefore, is the sum of jumps of the scattering phase. 



The same property takes place in almost all one-dimensional systems 
with any interaction. The reason is that in one-dimension any small bare 
interaction is effectively strong. Below we discuss the possibility that the 
topological order itself may be developed as a result of an effectively strong 
interaction in any spatial dimension. 

Several peculiar properties follow immediately from the assumption of 
a topological order. For example, one-particle excitations in the one- and 
three-dimensional quantum antiferromagnet have fermionic statistics, while 
in two dimensions the statistics is fractional (half fermionic). 

We start from reviewing the one dimensional case were topological char- 
acter of the ground state is the most transparent. 

2 Frohlich Superconductivity — One Dimen- 
sion 

Let us now review Frohlich's ideas (see e.g.,Ref. [|1^, [1^). We start from the 
Frohlich model of an incommensurate electron-phonon system with 



where and bq are electron and phonon operators, respectively. 

In one-dimensional electron-phonon systems the Peierls instability causes 
a lattice displacement 




(12) 



(13) 



q x,q 




(14) 
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and a modulation of the electronic density p ~ cos(2A;i;'X + if), where A ~ 
exp(— const/fyf^) << Ep is the value of the phonon condensate at 2/cf and ip 
is its phase . 

Frohlich noticed that the periodic density fluctuations of electrons are 
fixed only relative to the lattice and can easily travel with some velocity such 
that p ~ cos{2kp{x — vt) + ip). This could be compensated by changing the 
phase according to p = —2kpv. Therefore the current is j = pv = v{Nkp/7i), 

jx = , (15) 

and from continuity 

N , 

P = P0 + TT^ , (16) 
zvr 

where N is the degeneracy of an electronic state (spin, for example). These 
are Frohlich's equations. 

Assuming that the charge density wave (CDW) is not pinned, i.e., [^d'f — 
Vodl)p = —Ere, we obtain the eq. (0) with A^^ = ^^o, where fo = is the 
Fermi velocity. 

Since there is a gap in the electronic spectrum, and the only gapless mode 
reflects the ability of the CDW to move, Frohlich concluded that his system 
is superconductive. In fact the CDW is an ideal conductor rather than a 
superconductor due to pinning mechanisms. We already mentioned that the 
Landau criterion does not work in one dimension - there is no space to flow 
around an obstacle. 

Let us now derive the Frohlich equations formally. The first step is to 
pick out the fast variables and keep only the slow variables. In this case 
the slow variables are associated with electrons in the vicinity of two Fermi 
points ±/cf, 

a{x) ~ e^^^^^i + e-^^-^Vi? , (17) 
and phonons with momentum close to 

Y.bqe"'"' ~ —6^(2'=^^+^) . (18) 

q 9 

In the continuum limit we then obtain the so-called linear cr-model 



,2 



lAI |A 



|2 



g'^uj^ g^ 



+ i){ib - |A|e'^5'^)^ , (19) 
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where D = 7^(9^ + and 70,71,75 are two-dimensional Dirac matrices, 
and ZJ is a characteristic friequency of phonons. The modulus of the phonon 
field does not fluctuate much and is determined by its mean field value A ~ 
exp(— const/(7^) << Ep. This is the Peierls instability - as a result of the 
interaction a gap A has opened at the Fermi level. Let us stress that position 
of the gap is always at the Fermi level, so that the spectrum strongly depends 
on the number of particles (filling factor). 

Consider now a linear response to an electric field (there is no magnetic 
field in one dimension), which brings us to the subject of the axial current 
anomaly. The point is that the continuum model (|T9|) possesses a local axial 
gauge symmetry, iIj{x) ip{x) expi'y^alx), in additional to the ordinary 
gauge symmetry, iIj{x) ip{x) expia'{x). This would mean conservation of 
the axial current, <9^J^ = Vp + fo~^<9tj = as well as the ordinary current, 
dfij^ = dtp + Vj = 0. As a result of this "symmetry" an external electric 
field may be gauged away which in it turn means that the system does not 
respond to an electric field. This of course not true, due the axial anomaly 
- in the quantum theory the axial current is not conserved. 

The axial current anomaly tells us that an external electric field gives rise 
to the non-conservation of the axial current = ip'^^'^^ip, 

d,3l = ^i?^" (20) 

where we have set f = e = c = 1 and cZ; — > 00. Using that in ID, = e^j^^j^, 
we get a one- dimensional version of the (§). 

e,ud,j. = ^E^"-' . (21) 
/vr 

In particular this equation means that the electronic system is compress- 
ible, even though the single particle spectrum has a gap. Compressibility is 
an inherit feature of the Peierls instability in incommensurate systems - we 
already stressed that the gap always opens at the Fermi level. Therefore, if 
one adds a particle to the system, it will not go to the upper band to occupy 
the lowest empty state. Instead it rearranges the period of the CDW, so as 
to create one more level in the lower band. The energy of the system does 
not change much by adding a particle, so the system is compressible. 

This phenomenon implements the so-called level crossing. The fact is that 
in the presence of a kink in the phase, 0(cxd) — 0(— cxd) = 27r, the electronic 
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spectrum remains unchanged, except that one level appears at the top of the 
lower band with an energy E = —A. When we add a particle to the system, 
it will therefore create a kink in the spatial configuration of and an extra 
level to be absorbed in the lower band. The anomaly equation (^) reflects 
this phenomenon. We derive this well-known result further on. 

Strictly speaking, we cannot distinguish between solid and liquid in one 
dimension - there is no transversal current. Nevertheless, there is a global 
version of the Meissner effect - a current in a closed wire must satisfy 

J jdx = (22) 

where ^^^^ is the magnetic flux within the wire. 

Combining ( pO| , pTD with the Frohlich equations (|T^|T^) and using the re- 
lation jn = —dL/dA'j^^ we obtain a bosonized version of the incommensurate 
CDW: 

/V 2 
= j-ilidn^^r - -E'^^'v) ■ (23) 



4-/]- V 2 ■ ■ • TP 

Let us now turn to the commensurate CDW and consider an electronic 
system close to a half-filled band. At exact half filling, the CDW is two- 
fold commensurate and so the vacuum is two-fold degenerate. A canonical 
tight-binding model is 

+ h.c) + (24) 

n 

where A„ „+i is a fluctuating hopping amplitude and i^A is a phonon energy. 
In the continuum limit the half-filled case is described by the same a- model 
(17) but with a real A: 

L = ii{it)-i^^A{x))ij-% (25) 

9 

In this case there is no soft translational mode since the CDW is commensu- 
rate and is pinnned by the lattice. The excitations are kinks of A (x) which 
connect two-fold degenerate mean field vacua: A ±Ao when x —* ±oo. 
In the presence of the kink the electronic spectrum remains approximately 
unchanged, except for appearance of the so-called zero mode, a state with 
a zero energy, located exactly in the middle of the gap. This mode is 1/2 
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degenerate ( the index theorem), i.e. it can accommodate one-half of a parti- 
cle. Fractional degeneracy means that if the particle has an isospin, only one 
particle with spin, say up, can occupy the zero mode. In coordinate space 
the wave function is located in the core of the kink. The axial anomaly in 
this case tells us that the density of extra states is equal to the density of 
zeros of A{x) (see, e.g.,|lT^) 

p{x) = |5(A(x))^ . (26) 

Suppose now we dope the system by adding additional particles S. The 
system will lower its energy by creating the number of solitons (zeros in 
A(a;)) which is necessary to absorb all dopants. 

Due to the interaction between solitons they form a periodic lattice. As 
a result of that, zero modes develop a narrow band in the middle of the gap 
with a width of the order A ~ Aq exp(— const(5p/p) where Sp is the doping 
density and p the density of the undoped system. This band absorbs all the 
dopants and is always completely full. 

As in Frohlich's case, these solitons have a translation mode due to their 
topological origin: a soliton lattice can slide along the atomic lattice without 
dissipation. Let Xi be the zeros of A{x). Then the density of extra particles 
(dopants) is Sp{x) = — Xi). Displacement of the positions of zeros 

around their mean field values, Xi = Xi + (p{xi, t)/27r, give rise to fluctuations 
of the density 6p{x) = p{x) — p. According to (^) they obey the same 
Frohlich equations (0,|1^), 

= ^(^i^ud^¥> ■ (27) 
The anomaly equations (0), with an extra p, apply againQ 

e,M, = ^E^-"' . (28) 

Each twist of (f adds one additional state in the middle of the gap. Therefore, 
adding He additional particles gives rise to the topological charge of Ue = Q 
where Q = / • AH of this is true in incommensurate cases when the 

^ In the lower density limit the factor p is hidden in the velocity vq. Here and later on 
we neglect its space dependence. 
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system, after doping, has infinitely degenerate classical vacua. If the doping 
is a rational number, say, p/g, then the number of degenerate vacua is finite, 
namely q. The CDW is pinned, generally by an exponentially small potential 
A2(Ao/£/)^-2cosg0. 

Let us list some important messages that follow from this picture: 

(i) An energy gap in the electronic spectrum eliminates elastic scattering. 

(ii) The system has infinitely many degenerate vacua. They are character- 
ized by the number of solitons, i.e., by a topological charge \Q). 

(iii) A topological configuration ip{x, t) such as t = 0) = and / dx^pdx = 
27iatt—>- oo, which transforms one vacuum \Q) into another \Q + 1), is 
a low energy excitation with a dispersion uj{k) = Vok. Then "the whole 
system, electrons and solitons, can move through the system without 
being disturbed" 0. 

(iv) This, however, is not sufficient in dimensions greater than one. To 
achieve superconductivity in higher dimensions, it must have all these 
properties and in addition he a liquid. 

3 Fermionic Number and Solitons 

3.1 Level Crossing 

Let us consider fermions in a static (vector or scalar) potential. When the 
potential changes adiabatically, the fermionic energy levels also shift. If the 
chemical potential /i lies in the gap, then for an adiabatic change of the po- 
tential the levels generally cannot cross the energy level E = fi and therefore 
all occupied levels remain below /x. However, there are certain potentials 
which create some unoccupied levels below fi or force some occupied levels 
to cross the level of the chemical potential. This phenomenon is known as 
level crossing. Unoccupied levels appear below (or occupied levels appear 
above) the chemical potential and separated from the continuous part of the 
spectrum are what is called zero modes. 

There is an important theorem (the index theorem) which states that 
potentials which give rise to crossing levels, must have a topological charge. 
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We refer to them below as solitons. In Sect. ^ we discussed solitons (kinks) 
in one-dimensional models. Below we review them and discuss solitons in 
some models in 2 and 3 dimensions which will illustrate the essential physics 
of the topological mechanism. 

3.2 Models and Solitons 

3.2.1 One Dimension — Kinks. 

The most general Hamiltonian of the Peierls model is 

H = ajdx + f3ni + i757r2 (29) 

where the Dirac matrices a^:,/?, 75 = a(3 may chosen as the Pauli matrices: 
«x = CTs, P = (71,75 = —ia2- We also assume that the modulus of the vector 
(tti, TT2) takes a fixed value at infinity. 

The soliton here is a field 7f(x) which forms a homotopy class 7ri(S'^). Its 
topological charge is 

With a given value of the Ivf | at infinity Q is an integer number. 

3.2.2 Two Dimensions — Vortices. 

In dimensions greater than one it is necessary to include a vector gauge field 
in addition to a scalar field. 
We consider two models: 

A. QED3 

H= J2 ^iM^d^ + «m) + . (31) 

(7=1,2 

A soliton for this model is a vertex or magnetic flux with the charge 

Q = -f^dx. (32) 

B. Another model is non-abelian (SU(2)) QCD3 coupled to a meson field 

H = a^{td^ + A^t) + (3Tn . (33) 
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Solitons in this model are a combination of the hedgehog of the meson 
field (an element of the class 7r2(5''^)) and a vortex of the gauge field. Its 
topolological charge is 

^ f, TvF In ■ Dn X Dn.cPx , 
y ^ |7r| 4 |7r|3 Mtt ^ ^ 

where F — diA2 — ^2^1 + Ai x A2 and Dn — dif + 2A x n is the covariant 
derivative. A soliton with a topological charge of 1 has the asymptotic form 

7r3 = (7ri,7r2)^^ ^ . (35) 

at infinity. A general form of a vortex with an arbitrary topological charge 
is given by the homotopy class 7r2(S'^). 

Let % — 7ff/|7f| be an element of this class. We shall set 

Ai = ^Tr^iTT = AiT , (36) 

a vortex with a charge 

g = / tr(7rci7rci7r) . (37) 

A unit of the topological charge is carried by a zero of vf. The density of the 
topological charge is given by the "magnetic field" projected onto n 

^vortex ^ ^^vortex _ l^abc^aQ^^bg^^c (gg) 

SO that B^^'^^'^dx is a total flux. 

3.2.3 Three Dimensions — Monopoles. 

As a model in three spatial dimension we consider the SU(2) gauge theory 
with chiral bosons: 

H = a^{idf, + Afj^r) + i'y^Trr + m . (39) 

The spectrum of this model is symmetric since the matrix 75 anticommutes 
with the Hamiltonian. 
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Solitons for this model are the Polyakov - 't Hooft monopoles. At m = 0, 
and at r — > oo the monopole configuration is 

7f(r) \e,J-± . (40) 

Similar to the two-dimensional case the soliton configuration is given by 
( PUD and is an element of (712(5'^)). Then the monopole charge is carried by 
zeros of vf [|15[. It has the form of (^) where the integral now goes over a 
boundary surrounded zeros of |7r|. The density of the monopole's charge is 

^ = ■^£ijktr{ndi7tdj7r)uk (41) 

where 

= dk\m\^\) (42) 

is the density of zeros of vf. 

The "magnetic field" projected onto vf is a gauge invariant, 

= 2^ijk^^jk = ^ ^ijkTT OjTT (JklT . [4:6) 

So q = ^i/ji?™""^, which is equivalent to the Dirac equation for the monopole 
charge q = -^diBf^°'^. 



Rotation of the vector vr implies a gauge transformation [|T5|. It is con- 
venient to describe a multimonopole system in the unitary gauge. If is a 
transformation which rotates vf to the third axis, tt = g'^r^g, then the gauge 
transformation 

TT gng'^ = rg (44) 
A: ^ 9A,g-' - d,gg-' = Afr^ (45) 
converts the non-abelian monopole to the Dirac monopole: 

= -tT{T,d,gg~') . (46) 

This gauge is singular. It transfers a topological charge from the vr-field 
to the gauge field. A general form of the topological charge can be written 
in a compact form if one includes the mass m into the 4- vector 7Ca = {m, vr): 

V — / ^ijk^abcdi T^TT J" ^ • [^') 
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3.3 Anomalous Calculus 



In this section we calculate the so-called "fermionic number" or the Atiah- 
Patodi-Singer invariant rj for the models discussed in the previous section. 
It is simply the number of new states appearing above (or below) a certain 
energy when one switches on the potential. For Dirac particles it is the 
number of states appearing inside the gap. These states are refered to as 
zero modes or midgap states. Since without a potential the spectrum of 
the Dirac operator is symmetric, this number also measures the difference 
between the number of states with positive and negative energy t] = iV+ — iV_ 
(spectral asymmetry). If iV = A^_|_ + is the total number of levels, then 
the number of negative levels would be = y ~ l'"^- "^^^ "regularized" 
number of negative levels i.e. the change of the negative levels due to the 
potential, AA^ = A^_^ = —\r], is the fermionic number. 



The index theorem (see e.g. states that the fermionic number of the 
Dirac operator, or the number of extra states induced by a sohton, equals 
the topological charge of the soliton, AA^ = Q. Moreover the density of the 
extra states is equal to the topological density g(x): 

p{x) = {ij\x)ij{x) - ^) = q{x) . (48) 

We illustrate this general fact for the models listed in this section. 
The fermionic number for the Dirac operator may be defined as 

AA^ = -iTr[sgnif] . (49) 

This divergent quantity has to be regularized. One of the standard ways of 
regularization is to replace it by 

'''' = -T. h^wh^ 

The integrand now converges. 

3.3.1 One Dimension — Kinks. 

Let us start from the Peierls model (pOj), 

H = ajd^ + pTCi + Z757r2 = 0-3, P = ai, 75 = -ia2 • (51) 



15 



Then the square of the Hamiltonian is 

= (id,Y + + Eabaad^TTb, a,b= 1,2 . (52) 

Let us now expand the integrand of (|50|) in dn. In the first non-vanishing 
order we have 



2n J-oo \{id)^ + 71^ + [idy + 71^ + 2^) 

This is the only order which contributes to the fermion number. The trace 
in spinor space gives 

The trace in momentum space is conveniently calculated in the plane wave 
basis, 



Finally the integrals over p and z give 

2^1 



AN = -— I dx e^,^^ , (56) 



the topological charge (|30D of kinks. The number of extra states induced by 
a soliton therefore equals the topological charge of the soliton. 

The model with vri = deserves special interest. This is commensurate 
Peierls model {112 = A). It is defined by 

H = ajd^ + Z75A . (57) 

The spectrum of the Hamiltonian is symmetric. This means that if there 
is a state xI^e with an energy E then there is always a state jSip with the 
energy —E, except for E ^ 0. It follows from the fact that the Hamiltonian 
anticommute with the matrix f3. 

Now the soliton is a kink with A(— 00) = — A(+oo) = Aq. Setting tti = 
in (0), we obtain (|i2| ) in the form 

Q = -\j 5iA)dAdx . (58) 
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In this case an extra state appears in the middle of the gap (zero mode). 
The kink does not respect periodic boundary conditions. As a result of this 
a zero mode state may accommodate only 1/2 of the particle. In a system 
with a periodic boundary conditions kinks may appear only together with 
antikinks, so the total number of states remains integer. 



3.3.2 Two Dimensions — Vortices. 

A. QED3. We now consider the model ( pTD 

H= J2 i'lMtd^, + a^) + Pm]^^ . (59) 



(7=1,2 



. The square of its Hamiltonian is 

= {idf, + a^f + m^- (3{dia2 - SsOi) (60) 

Let us expand the integrand in (|y) in terms of P{dia2 — d2ai). The first and 
only non-vanishing order gives: 



""{id^ + a^y + m^ + zy ' ^^^^ 
The trace leaves the integral over p and z: 

nJ-00 J J (27r)2 [(p^ + a^)2 + m2 + z2]2 ^ ^ 

As a result we obtain (compare (p2D) the relation between the fermion number 
and the topological charge of the vortices 

AA^ = — ^ sgn m [ d^x {dia2 — d2ai) . (63) 

B. QCD3 coupled with a meson field. This model is defined by 

n = a^{id^ + A^t) + Pnf (64) 
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and the square of the Hamiltonian is 

= {id^ + A^rf - (3Ft+ (7f)2 + a^p id^n r + a^(3i[A^ x 7r]2r . (65) 

To simphfy calculations we first drop the gauge field. Then, expanding 
the integrand in ( |50| ) in terms of id^n f up to the second term, we have 



1 /• + 00 

AA^ = / rf^ Tr 

2tt 



CO 



1 1 - {ttf^pid^n f) ^ — - + ( {a^pidf^n f) 



id^Y + (7r)2 + ^2 ' M J ^^Q^^2 + (^)2 + ^2 



(66) 



Changing to and taking the trace in the spinor and isospinor spaces 
we get 

^^ = -2^i-oo'^i'" y(2^ [(p,)2 + (vf)2 + .2]3 ^/^[^^^^- ^2vrr] . 

(67) 

Finally, integrating over p and z we obtain the fermion number 



AN = jp{x)dPx (68) 

with the density 

P = -e--e-" ""^;!'f-"^ . (69) 

87r |7r|3 ^ ^ 

Once again we obtain that the density of the extra states is equal to the 
topological charge of the soliton. 

To include the gauge field one must first replace derivatives in (|69|) by 



covariant derivatives. Then we notice that for tti = = and Aj^ = A'j^ = 
the fermion number must reproduce the QED3 result (|63|) . Taken together 
these considerations give 

^ 1 ^^.^atc ^^d.n" + 2[A X nf)id,n^ + 2[A, x tc^) 

^ 2tt \tt\ 8tt |7f|3 ^ ^ 

where F = 81X2 - 82^1 + [Ai x A2]. 



18 



A special case of interest is tts = 0,Ajj^ = Aj^^ = 0. In this case the 
spectrum is even as the matrix /3r3 anticommutes with the Hamiltonian. 
Equation ( [7D| ) gives 

p = -ieijdi7T,dj7T2 + i^|v?|'V X A'')5{n) . (71) 
47r 

In this case all extra states are located exactly in the middle of the gap (zero 
modes). Each zero mode may accommodate 1/2 of a particle. 



3.3.3 Three Dimensions - Monopoles. 



The fermion number for the model ( PU] ) is again equal to the topological 
charge of the monopole { ^7\) ||16| . The simplest way to calculate it is to drop 



the gauge field and proceed similarly as in the previous section. Now we 
must keep the third order term in the expansion ([66|). As a result we obtain 

1 /■ ,3 TTadiTCbdjTfcdkTCd 
2lT^ J TT 



= — / d^xeijkEabcd pn • (72) 



Next we want to make this expression gauge invariant. After a replacement 
of derivatives by covariant derivatives, it still remains gauge non- invariant. 
To rectify this one must add the second term in (^7]). 

Similar to the previous examples a very special case arises when m = 
ttq = 0. In this case the matrix (3^^ anticommutes with the Hamiltonian, 
and the spectrum is symmetrical. Therefore the extra states appear in the 
middle of the gap at E = (zero modes). As usual each zero mode is 1/2 
degenerate. Setting ttq = in (^7]) we obtain 



1 1 

P(^) = ^£ijktT{Tf[di7cdjfc + —Fij])uk (73) 
ovr 4:71 

where it = 7f/|7f| and Uk = dk\TT\6{\Tf\) is the density of zeros of vf. 



4 Topological Fluid as a Superfluid 

4.1 The Basic Model — Holons and Spinons — Zero Cur- 
rent Theory 

The main inspiration for considering topological superconductivity comes 
from strongly correlated electronic systems. They suggest a general model 
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to study this phenomenon. Running ahead we discuss this model now and 
later derive it from a canonical model of the doped Mott insulator. 

Consider two sorts of particles "holons" h and "spinous" (the names 
come from the Mott insulator literature) which obey the so called "zero cur- 
rent" theory. This means that densities and currents obey the local constraint 

(i) their densities are always complementary, h^{x)h{x) = il!{x)ip^{x), namely 
the number of holes in spinous is the number of holons, where the total 
number of holons is given by / h'^ {x)h{x)dx = 6 (doping) and 

(ii) their currents are always opposite (the so-called "zero current theory" ) 
Jhix) +jsix) = 0. 

The holons carry an electric charge and spinous are neutral, so that the holon 
current is in fact the electromagnetic current. The Hamiltonian of the basic 
model is 

n = —hUx)H^h(x) + tlj(x)H^(x) + , (74) 

m 

where H may be any Hamiltonian which exhibits a level crossing. For ex- 
ample, one may consider any one of the models considered in the previous 
section. The local constraints are automatically enforced if H contains an 
abelian gauge field, as in (^) of QED3 , other wise a gauge field ( Lagrangian 
multiplier) must be added to (ff^): A{jh + js) + Ao^h"^ {x)h{x) — 4'{x)ip'^{x)). 
Later we derive the Hamiltonian H for the Mott insulator. This Hamiltonian 
is in fact different from any of the standart models listed in the Sect. |^, but 
has essential common features with models of the sec. 3. 

4.2 Topological Instability 

Let us consider the basic model of topological fluids (0). We argue that the 
perturbative vacuum of the model is unstable with respect to the creation 
of topological charge and solitons, namely we show that the vacuum has 
a topological charge such that the fermionic number equals the number of 
dopants (holons), 

-^7] = AN = S . (75) 
This the topological charge of the ground state. 
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To see this, let us apply the adiabatic strategy: first consider all poten- 
tials A and vr static and find the best spatial configuration to minimize the 
energy for a given doping, and then consider fiuctuations around the minimal 
static configuration. First assume that there are no solitons and start adding 
dopants. Consider first the spinon part of the model. Dopants start to fill 
the upper band of the Dirac spectrum of spinous. This costs the energy of 
the gap, plus the Fermi energy of spinous, plus some small radiative correc- 
tions of order of 6. If, however, there are enough solitons to absorb by their 
zero modes all dopants, we pay less than the gap and gain all that energy. 
In addition we gain the Fermi energy, since the zero modes form a narrow 
band inside the gap. The same arguments hold for the holon part of the 
Hamiltonian as well. The solitons form a narrow band for the holons at zero 
energy, so they also favour the topological vacuum. As a result the number 
of solitons in the vacuum must provide that number of zero modes which will 
absorb all dopants, h'^h = 6. 

To formalize these arguments, one can show that the perturbative vac- 
uum is indeed unstable. The topological condensate appears in the one-loop 
correction. In leading order of the chemical potential, the free energy is 

E -l/2\fxri\- fi5 (76) 

where E is a perturbative energy. Equations ( |75D minimize this free energy. 

We conclude that a small doping requires a configuration of vr and A with 
the relevant topological charge determined by (fTSl). In turn the wave function 
of the state with a topological charge establishes the Chern class (7) equal 
to the charge. We refer to this phenomenon as topological instability. 

4.3 Topological Field Theory coupled with the Matter 
Field and Superfluidity 

A simple version of the basic theory occurs when the gap in the Dirac spec- 
trum is the largest parameter in the system. Integrating over spinous we can 
then neglect all terms except the so called the Wess-Zumino-Witten term, 
denoted WZW hereafter. It is the only term which does not vanish when the 
gap approaches infinity, and is given by 

WZW= lim InBetiidt + H) (77) 



21 



The WZW term captures in the Lagrangian language the phenomenon of 
anomahes, namely, anomalous currents may be obtained as the Lagrangian 
equation = — (5WZW/(5yl^. For the ID model (pUD, the topological current 
(see dO)) 

Therefore, 

WZW = - / pe^^SatA,^^^ (79) 
, The WZW term for QED3 it is known as the Chern-Simons term: 

— sgn me,, xA^d, Ax . (80) 

For QCD3 and for the 3D model WZW = — / A^q^dx, where is the topo- 
logical current given by the integrand of (32) and (42) |[T7[] . 

The WZW term itself gives rise to and is descriptive of topological field 
theories. A topological field theory describes only the zero mode sector of a 
field theory and the zero mode sector is invariant under all diffeomorphisms of 
the volume - positions of the solitons are not fixed unless they interact with 
matter. This means that any external perturbation, such as inhomogeneous 
pressure or twist, does not change the energy of the system. As a result the 
zero mode sector is topological in the sense that it has no physical states in 
the bulk. Physical states appear only as a result of a boundary and they are 
confined to a boundary (the edge states). 

Our basic model is a topological field theory coupled with matter (holons), 

L = J dx h\idt + Ao + H^)h + WZW . (81) 

In this case the zero modes are occupied by particles (holons) and possess 
some dynamics. Nevertheless, this dynamics is limited as a result of the 
overall gauge symmetry of the unperturbed topological field theory. Matter 
of course destroys all diffeomorphisms which change the density of the matter, 
so that only diffeomorphisms which preserve the density will survive - but 
matter which is invariant under density preserving diffeomorphisms is an 
ideal liquid i.e. a superfluid. 
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4.4 Superconductivity 

As in the one-dimensional case, it is almost obvious that a dilute incom- 
mensurate system of solitons forms a compressible liquid. Indeed, a soliton 
position is not fixed relatively to the rest of the system. It can be translated 
to any place in the system without changing the energy - zero modes are 
translationally degenerate. Therefore, if a soliton moves slowly, say with a 
momentum k, the energy of the system changes by vk. Since the interaction 
between solitons is short range, the same remains true for dilute system of 
solitons Q Thus we assume the solitons form a liquid. Each soliton carries a 
density of charge carriers and an electric charge according to its topological 
charge ([75|). Translation degeneracy of the zero modes therefore implies lon- 
gitudinal sound in modulation of electronic density and this conforms to the 
first part of the Landau criterion. 

Furthermore, the narrow band formed by electrons trapped in solitons 
is always filled and separated from the rest of the spectrum. Therefore all 
excitations, except coherent modulation of charge density and topological 
density, cost non-zero energy. Indeed, one particle excitations consist of a 
particle being taken away from the soliton core. Since all particles form 
bound states with solitons, this costs energy i.e. the excitation has a gap. 
This conforms to the second part of the Landau criterion. To summarize, we 
conclude that the Landau criterion is satisfied and the topological liquid is a 
superfluid. 

To formalize these arguments let us consider for example the basic model 
(|75) with 

H= J2 + «m) + (82) 

o-=l,2 

in two dimensions. At low energy the system develops a density of magnetic 
flux determined by the density of extra particles: 

= --V X a = --F(x) . (83) 

TT 71 

^ At larger density solitons may form a Wigner crystal, i.e. a solid, or the interaction 
may destroy a bound state between a soliton and electrons (whatever comes first). It may 
happen that the soliton density is commensurate with the lattice. In this case solitons 
form a solid and arc pinned by the lattice. If commcnsurability is weak, the solid may be 
melted by quantum fluctuations. We do not discuss these questions here. 
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When the sohtons move, they are followed by particles, implying 

f[x) = --{dta-Va^) = --E , (84) 

where j* = Eikjk, and an additional factor 2 comes from the spin. A 
homogeneous density p gives rise to a homogeneous magnetic field (F) and 
no electric field. This is the first level of adiabatic approximation. To study 
modulation of the density and to see that there is no dissipation, we want 
to see how the gauge field fluctuates. Its fluctuations are determined by the 
motion of electrons in the presence of magnetic flux F/ir = p. The relation 
between the density and topology implies that the first two Landau levels are 
completely filled and in addition they separated by the gap m from the upper 
Dirac band. In other words, all single particle excitations have a gap, which 
means that excitations over the gap produce an energy ^eE^ + -^{F — {F))"^. 
Here e and x are the dielectric constant and diamagnetic susceptibility of 
electronic zero modes, respectively. As a result static fluctuations of the 
topological charge F{x) — (F) and topological current E{x) are short range: 

{{F{x) - (F)), (F(x') - (F))) = xSix - x') , (85) 
{E{x), E{x')) = e-^6{x - x') . (86) 

We can now obtain density and current correlation functions since they 
are directly related to fluctuations of the topological charge. Combining them 
with the continuity condition, we obtain the linear response function for a 
superfluid with = x and v"^ = ex in (P and @. 

A similar strategy applies for the 2D model (|33D . Analysis of the three 
dimensional model is more complicated and involves some new features [^]. 
The physical picture, however, remains the same, as briefly indicated below. 

Since all single particle excitations have a gap, we expect that all low en- 
ergy excitations are fluctuations of the topological charge go and topological 
current g^. For the same reason, we know that static fluctuations are short 
range. We also know that there is a longitudinal sound since the topological 
configuration possesses a soft translational mode. Assuming that solitons 
do not form a crystal, we must conclude that the low energy topological 
excitations obey hydrodynamics: 

H = fyo+Vo'f] ■ (87) 
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where {qo,<l) is a density of the topological charge and topological current of 
solitons, 

Finally, as a result of the topological instability, the density of electric 
charge (current) is directly related to the topological charge (current) 

= '^It^ ' (88) 

which, together with (0), gives the hydrodynamics of the ideal charged 
liquid, 

. i.e. hydrodjTiamics of the superconductor (|^). 

4.5 Topological Superconductivity and Bosonization in 
Higher Dimensions 

Bosonization is a fascinating phenomenon of one-dimensional physics - col- 
lective modes of interactive fermionic systems may be described by bosonic 
particles. Moreover, fermions themselves can be treated as solitons of bosonic 
non-linear waves. The physical reasons behind bosonization are the same as 
for the Frohlich superconductivity. In fact, a generic one-dimensional incom- 
mensurate system develops some density wave (CDW,SDW, etc.) as a low 
energy excitation. These are bosonized fermionic collective modes with 

H = + vl{dM'] ■ (89) 

(Here vr is the canonical momentum of the bosonic field ip.) In other words, 
bosonization in one dimension implies hydrodynamics. If the one-dimensional 
system is commensurate, the Hamiltonian of the density wave acquires some 
non-linear terms, as in the Sine-Gordon model, for example. 

This interpretation allows the extension of bosonization into higher di- 
mensions. We shall say that a fermionic system can be bosonized if its low 
energy excitations can be described entirely in terms of currents, such as 
when an effective Hamiltonian of low energy dynamics has the Sugawara 
form, i.e. may be written as a bilinear form of currents. Recall that this 
form implies the hydrodynamics of the superfluid (6). The Sugawara form is 
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a linear form of quantum hydrodynamics [|^ : 




J = pv . 



(90) 



The bosonized version of hnear hydrodynamics appears in terms of a 
displacement u. Let us solve the continuity condition by setting p — p = pS/u 
and i = pdtu. Then (6) becomes 



In contrast to one dimension, not every two-or three-dimensional fermionic 
system can be bosonized, but only those which are superfiuids. For instance 
a fundamental property of the Fermi liquid is that it cannot be bosonized. 
Dissipation and gapless excitons are its own characteristic features. 

We have indicated that a fermionic theory with a topological charge in 
the ground state can be bosonized, regardless of its dimension. From this 
point of view bosonization can be treated as a projection onto a zero mode 
narrow band. 

5 A Gauge Model for Correlated Electronic 
Systems 

Strongly correlated electronic systems represent physical systems where topo- 
logical fluids may appear. In this section we sketch the derivation of a basic 
model of a topological fluid from the Hubbard model or the "t - J" model - 
a canonical model of strongly correlated electronic systems. We also empha- 
size a set of assumptions and approximations which are necessary to separate 
the physics of topological fluids from a variety of other physical phenomena 
which occur in a correlated electronic system. 

5.1 A Model for the Doped Mott Insulator — Resonant 
Valence Bonds and Chirality 

The "t — J" model is defined by 




(91) 



H — '^^^ij^ta^jc! "I" JijSiSj (all Ji j > 0) , 



(92) 
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where Si = Cia(Taa'Cia' is a spin operator of an electron at the lattice site i. 
A hopping amplitude tij and an exchange amplitude Jij connect the nearest 
sites and the total number of electrons is close to the number of lattice sites: 
Ne = Nq{1 — p), while a strong Coulomb interaction does not allow doubly 
occupied states: 

n^= J2 c+c,^ = or 1 . (93) 
'^=1,1 

At zero doping the model describes the Heisenberg antiferromagnet 

H = J2JabSaSt , (94) 
ab 

with spins Sa and Sb on sublattices A and B, respectively. 

In addition to an average value of spin (classical description) there are 
two operators which characterize the ground state of an antiferromagnet, 
namely density of energy 

e,, = (1 + Sjj) (95) 
and chirality or measure of topological order jn], ^ Q 

WiC)=til[il/2 + aS,) , (96) 

ieC 

where a are Pauli matrices, and C is a lattice contour. 

The latter operator is of particular importantance for the doped case since 
it determines the correlation of electronic phases at different spatial points. 

These two operators are related to the amplitude and phase of An- 
derson's Resonance Valence Bond (RVB) through 

e^, = \A,,\' (97) 

and 

WiC)=l[A,, . (98) 

c 

It follows from this definition that the RVB is a gauge field. It can be locally 
transformed by a U(l) transformation 

A,, ^ A,,e^(--"^) (99) 



27 



and by a symmetry group of the crystal class 

A ^ C-^/\C , (100) 

without changing the physical properties. In terms of this field, the topolog- 
ical order parameter W{C) acquires the form of a lattice Wilson loop: 

W{C) = e**(^) . (101) 

Roughly speaking, the flux of the RVB field 



e 



Y[e'^^J , (102) 



c 



where Aij is a phase of Aij, represents a "magnetic" flux, and penetrates 
through a surface enclosed by the contour C. The RVB field can be formally 
introduced by fermionic representation of a spin operator: 



with the condition 



so that 



S = cffffaa'Cia' (103) 

Ec^tci. = l (104) 



A., =Ec.W . (105) 
Using this representation the Heisenberg Hamiltonian can be explicitly 



written in terms of the RVB field 18, 19 



as 



Ha= 4^U^,c„J + ^o{^){ct,c„^ " 1) + U,,JJ^^Uh ■ (106) 

(T, (ab) 



The field Aq is a Lagrange multiplier which ensures the constraint (p3|), while 
the phase of the conjugated RVB field Uij forms a time and space component 
of a gauge field. The constraint ( |104|) is relaxed after doping, when a small 



density of empty sites appear on the lattice, 

Y.4^c„. = ^-hthi , (107) 
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where hi is a dopant (hole) operator. An electron operator is then represented 
by the product of a spinon and a holon, C(jih\. The hopping Hamiltonian is 
therefore 

H ^tJ2hHa)h{b)A^^^+ c.c. (108) 

{a,b) 

where a and b are sites of sublattices A and B and A- ^ = {J2a ct(^)co-(^))- 
So far all manipulations were exact. 

Depending on parameters this model may describe physically different 
situations. We are interested in a situation when fluctuations of the bond 
energy (modulus of Aj ,,) are much less and much slower that fluctuations of 
the bond phase (chirahty ) . If their scales are separated we refer to this state 
as quantum antiferromagnetism. Let us stress that in the Necl state the mod- 
ulus and phase of bonds fluctuate similarly and cannot be separated. In this 
case holes form the Fermi-liquid and interact weakly with spin fluctuations. 



5.1.1 Adiabatic Approximation. 

Some progress can be made if the magnetic background can be treated adi- 
abatically. This means that a hole moves in a slowly varying spin configu- 
ration, similar to the Peierls-Frohlich problem where electrons move with a 
slowly varying distortion. In this case we may use a semi-classical strategy: 
first find a static A and Aq which minimize the energy of the system with a 
given doping, 

and then take into account quantum fluctuations around the static mean 
field. 

To apply the semi-classical strategy one more step is required. The point 
is that even if the hopping amplitude is bigger than the exchange energy, the 
dynamics of spins cannot be considered adiabatically - in the antiferromagnet 
each jump of a hole sharply changes the spin configuration by fiipping a spin 
on a sublattice. However, two consecutive jumps bring a hole to the same sub- 
lattice, so the spin configuration remains approximately unchanged. Thus, 
to apply adiabatic arguments we must first integrate over fast and sharp pro- 
cesses. Perhaps, the easiest way to implement the adiabatic approximation 
is to introduce a difference between energies of a hole on different sublattices 
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by adding the term fJ^iJ^dCgCg — Z^gclc^) to the Hamihonian. If the hopping 
energy t is less than fi ( adiabatic parameter), a holon and a doublon appear 
on different sublattices only virtually. Therefore in leading order in t //i one 
may consider only processes of two consecutive hoppings described by the 
effective holon Hamiltonian 

H = t'{Y: h\a)^l,,h{a') + h\h)^l^,h{h') + c.c.) (110) 



where 



= E \b'^h' = E \d'K,y (111) 



are hopping operators within the same sublattice. After all, we suggest that 
this Hamiltonian be considered as phenomenological. 

Now as long as t' ~ >> J, one can treat magnetic fluctuations 

adiabatically. A solution of the mean field equation ( |109D depends on a 



particular exchange coupling. Here we are interested in a solution where 
|Aafe| slightly varies from bond to bond, but never vanishes: A^b 7^ 0. Then 
the mean field values of Ajj and Uij are proportional to each other. After all 
these considerations we end up with the basic Hamiltonian 

b,b' a,b 

+ J2Ao{clc^^ + hlh,-l) , (112) 

i 

where t ~ t'J~^ and J ^ J are phenomenological constants and we set 
|Ajj| = 1. To remain within adiabatic approximation we also assume that 
i> J. 

We refer this model as the gauge model of correlated electronic systems. 
In addition to holons and spinous it has a gauge field, the phase of the bond 
Ay = expiAij, and a Lagrange multiplier Aq as a dynamical gauge field. Let 
us also notice that the gauge field has no energy of its own. 

The hopping part of the model has an important symmetry which reflects 
the symmetry of the square lattice. The spectrum of the Hamiltonian is 
symmetric, namely, if there is an eigenstate with an energy E, then there 
is an eigenstate with the energy —E. To see this one must note that the 
Hamiltonian ( [L12| ) changes sign if one replaces cg and /ig to — cg and — /ig and 
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keeps and unchanged. It also means that there is a matrix T which 
anticommutes with the Hamiltonian. In the sublattice basis it is a diagonal 
matrix 

r = diag(l,-l) . (113) 



5.2 Flux Phase 

The next step in the adiabatic approximation is to determine a mean value 
of the gauge field Ay, Aq. We approach the problem in two stages. At first, 
we neglect a small doping and find the most favorable configuration of Uij 
for the half-filled case. Then we discuss the effect of doping. The first stage 
is again approached in two parts. Assuming that \Aab\ is a constant, we first 
find the phase of Aab or Uab i.e., the ground state value of the chirality (^) 
and later consider a small variation of \Aab\- 

At the first step, we suppose that a variation of \Aab\ in space is small 
and may be found on the basis of continuous theory. Let us therefore set 
\Uab\ = AoJab = tg in the first approximation. Then the only degree of 
freedom to be determined is the chirality W{C) = Aq^ ' exp{i(f){C)). It has 



been suggested in Refs. E^, pOl 21| and has been proved recently by E.Lieb 



that the fiux vr per elementary plaquette provides the minimum of the 
energy: 

WiP) = A^(-l) . (114) 
This state is known clS db fiux phase or chiral phase. 



5.3 Continuum Limit 

As usual, in order to go to the continuum limit we must separate fast and 
slow variables. The fast variables are represented by the gauge potential 
of the background fiux. In the fiux state translations do not commute any 
longer. They form a group of magnetic translations: 

T,T,^+T,T,^ = , (115) 

where Cj are primitive lattice vectors. Apart from spinor representations of 
SO3, it is also reprented by 2'^ matrices of dimension d x d (where d is a. 
spatial dimension) which obey the Clifford algebra 

{ai,aj} = 26ij . (116) 
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Setting Aab to its mean field value Aab = a^, and Ao{a) = Ao(6) = m, 
we get the nearest neighbour part of the Hamiltonian in the form 

d ^ d 

H = t^T^^e^^^*^^"') + mp = t^aiCos2fci + mp (117) 

1=1 i=l 

where the wave vector k = {k^, ky, ...) belongs to the reduced Brillouin zone 
— I < fcj < ^, while the matrix /3 = ±1 depending on whether a site belongs 
to sublattice A or B, and = 0. To find the spectrum it is sufficient 

to consider the square of H. Due to (|115D one has 



= Y,cos2ki^ + . (118) 

1=1 

The spectrum contains two Dirac bands which touch at kp = (|, |, ...). 
The lower band is completely filled. Expanding in small momentum p = 
k — kp around the Dirac point, we get the Dirac Hamiltonian 

H = Y,i^l{c^d + m) (119) 

where the Dirac spinor ip^j is the slowly varying part of Cg-. To see this, let 
us denote a lattice site R as 

3 

R = Y.{2x, + (r, - l))ei , (120) 

i=l 

where rj = 1, 2 and Xj are integers. If rj = 1, the site belongs to the sublattice 
A] if Tj = 2 then to a sublattice B. Then 

i^vA^) ~ exp{ik„R)cj^^^ . (121) 

Let us mention that the dimension of the matrices a is twice that of 
the irreducible representation of the Clifford algebra, so they belong to a 
reducible representation. This is the well known lattice doubling - a lattice 
always produces an even number of fermionic species. To summarize, let us 
denote Tg^ = af, 70 = ^; 7, = Then {a^, aj} = 25^; (3ai (3 = af, = 
1 and the 7^'s form an algebra of Dirac matrices, {7^ 1u} = '^dfiu ■ 

This type of mean field theory is essentially invalid, because it does not 
manifestly reflect the main physical property of the theory - its gauge invari- 
ance. However, the gauge symmetry is restored by fluctuations. Before we 
discuss gauge fluctuations around the mean field, let us choose a convenient 
basis of the a matrices on a square lattice. 



32 



5.4 Lattice adopted basis 

Expansion around the mean field value follows the strategy developed in 
polyacetylene physics |jl2|, |l3l - only processes with a momentum transfer of 



2kp are relevant. Others lead to terms with higher order derivatives. For the 
sake of simplicity we consider here the basis adopted on the 2D square at half 
filling. Generalization to the three-dimensional cubic lattice is developed in 
Ref. @]. Since at half filling the period is doubled, we must consider the four 
nearest points of a lattice sell as a set of discrete variables. Let us denote 
them in terms of Cartesian coordinates: r = {r^jVy) : rj = 1, 2. Then a can 
be chosen in the form 

ay = l(g)ai; /3 = (Ti ® 1 , (123) 

where the cr's are Pauli matrices. In addition there are three more matrices 
Tx, Ty and Tz which commute with the matrices a. We choose them to form 
a basis of SU(2), [t\ t^] = ieijkr'': 

Ti = ail ; r2 = 0-3 (g) o"! ; rs = (72 (g) (Ji . (124) 

There is of course a direct product basis in which spin and isospin appear as 
in independent spaces: 

a ^ MaM"^ = 1 O ai ; P 10 p ; r^rOl. (125) 

In continuous theory, the amplitude Aab fluctuates around its mean field 
value such that these fluctuations are smooth, but not necessarily small. 
Fluctuations of the modulus of Aab are small and we neglect them, but 
fluctuations of the phase of Aab restore the gauge invariance and may not be 
small. To describe them we write 

Aab = e^^aiTi , (126) 

where Tj = exp{—di) is a translation operator and Aj a diagonal matrix. 
(Fluctuations may change the value of a bond, but do not introduce a new 
bond.) Four diagonal matrices which form a basis for fluctuations may be 
written in terms of the Dirac basis. They are 1 ® H and given by 

(3t^ = a^^ffS _ (^27) 
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The simplest way to include fluctuations of the phase is to restore the 
gauge invariance. The original lattice model is invariant under the U{1) 
gauge transformation Cj Ciexpiipi. In the continuum limit where each of 
the four sites is treated separately, the gauge group becomes a four parameter 
group, and accordingly the phase of the gauge transformation is an arbitrary 
diagonal matrix: 

ip e''^^^^''i) . (128) 

To restore the gauge invariance of the mean field Hamiltonian ( |119|) we 
must add compensating fields: 

H = exp[ii?i/9r^ + iGiaiT^ + iKia2T'^]ai{idi + ai) exp[— ■iE'i/^r^] 
+ exp[i£^2/3r^ + iG2aiT^ + iK2a2r'^]a2{id2 + 02) exp[-z£'2/5r^] 
+ pT^m . (129) 



These compensating fields are gauge fields. Under the transformation ( |12^ ) 
they transforms as follows: 

G^ + ip2 
+ LPs 

^ E^ + if^ (130) 

The quantities E = Ei - E2,G = d - G2, K = Ki - K2 and F = 
V X a are therefore gauge invariant objects. They are not small, so the 
exponents in ( |1 29|) cannot be expanded. The four types of "magnetic field" 



describe fluctuations of the flux in four adjacent plaquettes and implement 
the symmetry of the square lattice. 

Fixing a gauge, one may slightly simplify the Hamiltonian (|129|) . One 



may choose G2 = 0, -ft'i = 0, ii^2 = 0, in which case the Hamiltonian has the 
form 

H = exp[iE f3T^]ai{idi + ai) exp[—iE (3T^]a2{id2 + 02) exp[—iE2(3T^] 

+diGT^ + d2KT'^ + pT^m (131) 
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We note that even though the terms idiG and appear to be temporal 
components of the SU(2) gauge field, their symmetry is different. They are 
components of the spinor connection of the original crystal group. 

Fluctuations of the modulus can also be obtained from this symmetry. 
First of all the Hamiltonian must anticommute with the matrix (|113|) , which 
in our basis is F = (3t^. Secondly, since the Hamiltonian connects only the 
nearest lattice sites, it must be a tensor product of diagonal and anti-diagonal 
matrices in the sublattice basis. Except for the matrices a^, OiT^, r^, T2 which 
already appeared in the Hamiltonian, three more such matrices, /3ajr^, f3T^, ^t"^, 
describe fluctuations of the modulus of Aab- Taking them into account too, 
flnally gives 

H = exp[iE j3T^]ai{idi + ai) exp[—iE PT^]a2{id2 + 

+diGT^ + d2KT^ + i(3aiT^^i + /5r'*<l>'* + (3T^m (132) 

where $j, $° are components of the modulus of A in a proper basis. As a 
result, the basic model for the doped Mott insulator in the flux phase has a 
form: 

n = —hUx)H^h(x) + ^Ux)Hij(x) + ^($2 + (<I>^Y) . (133) 
2m 2A 

For further explicit reference, let us present the Hamiltonian H for the 
three dimensional flux phase. In this case translations are 8x8 matrices. 
They can be chosen in the form 

Ti = tti = (Ti ® (T3 (g) 1 

T2 = a2 = 1 (S) cTi (S) cTs 

T3 = as = as O 1 (g) 0-1 . (134) 

The fourth Dirac matrix f3 may be chosen to be diagonal, 

= (X3 ® (T3 ® (T3 . (135) 

Four other matrices which commute with the Dirac matrices are translations 
along face diagonals, 

T12 = = (72 ® o"! ® 0-3 

^23 = = (T3 (g) 0-2 ® (Ti 

T31 = = C7i ® (T3 ® (Ta , (136) 
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and translation along the space diagonal, 

Ti23 = hblo = 7l7273 = «ia2tt3/5 



0"! ® 0"i ® 0"i 



(137) 



This last operator is known as the parity operator. 

The 8 diagonal matrices which describe fluctuations are 



Oi = 1 ® 1 ® 1 

O3 = era ® 1 ® 1 
O5 = 1 ® 1 ® 

O7 = 0-3 (g) 1 ® 0-3 



O2 = 0-3 ® 0-3 (T3 

O4 = 1 ® (Tg (g) 1 

Og = 0-3 (g 0-3 1 
Os = 1 ® 0-3 ® 0-3 • 



:i38) 



Following the procedure we discussed for the 2D case, one obtains 
H 



^exp \ iEip + {iGiT^amanP + iK^T^aynan)\oii{idi + a^) x 

X exp iEi(3 + ^ (iG^T^amCHnP + iK^r^amdn 



(139) 



Despite some efforts, topology of the numerous fields appeared in the 
models ( [133|jl39| ) and their anomalous properties have not been elaborated. 
Nevertheless, there is some confidence that it will exhibit topological insta- 
bility and topological superconductivity along the scenario outlined in these 
lectures. 



6 Conclusion 

In summary, we have seen that if the ground state of a quantum system 
dynamically acquires a non-zero Chern's number (i.e. a topological order), 
then it is superconductive. In units where e = c = m = 1, the London depth 
is literally the density of the topological charge. This connection is rigid and 
as is the case with the whole theory, it is determined by the geometry of the 
phase space. It therefore seems likely that superfluidity can be derived on a 
more phenomenological (model independent) basis. 

It is also clear that topological superconductivity is drastically different 
from the BCS mechanism. How does this difference manifest itself in ob- 
servables? For example, what are the tunneling properties of the topological 
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superconductor, and the symmetry of the order parameter? This is perhaps 
the most interesting avenue of developing this theory. 
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